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Summary. — A report on the recent next-to-leading order QCD calculations to 
tibb and tijj at the CERN Large Hardon Collider is given. The elements of the 
calculation are briefly summarized and results for integrated and differential cross 
sections are presented. 
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1. Introduction 

Even though first results for next-to-leadiiig order (NLO) QCD corrections to heavy 
quark production were presented in the late 80's and early 90's [1-4] the topic of higher 
order corrections to tt is still very active and far from complete. Our present level 
of understanding is very well summarized in experimental and theoretical reviews, see 
e.g. [5-7]. Recent progress in NLO [8,9] and next-to-next-to leading order (NNLO) 
[10-18] calculations, as well as next-to-next-to-leading-log resummations (NNLL) [19-21] 
for inclusive tt hadroproduction is truly astonishing. The list for the more exclusive 
channels is just as impressive: NLO QCD corrections have been calculated for the ttH 
signal [22-27], where the Higgs boson has been treated as a stable particle. Most recently 
the factorisable QCD corrections to this process have been presented [28], where higher 
order corrections to both production and decay of the Higgs boson into a bb pair have 
been calculated. Moreover, NLO QCD corrections to a variety of 2 ~> 3 backgrounds 
processes tij [29-31], ttZ [32] and ttj [33] have been obtained. Last but not least, the 
NLO QCD corrections to the 2 — >■ 4 backgrounds processes such as ttbb [34-37] and 
ttjj [38] have also recently been completed. 

Both processes pp — >■ tibb and pp — > ttjj represent very important background reac- 
tions to searches at the LHC, in particular to ttH production, where the Higgs boson 
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decays into a bb pair. A successful analysis of this particular production channel re- 
quires the knowledge of direct ttbb and ttjj production at NLO in QCD [39]. In this 
contribution, a brief report on these computations is given. 

2. Theoretical framework 

NLO QCD corrections have been calculated within the Helac-Nlo framework. It 
consists of Helac-Phegas [40-42], which has, on its own, already been extensively 
used and tested in phenomenological studies see e.g [43-46]. Helac-Phegas is a multi- 
purpose, tree-level event generator which is the only existing implementation of the 
algorithm based on Dyson-Schwinger equations. It can be used to efficiently obtain 
helicity amplitudes and total cross sections for arbitrary multiparticle processes in the 
Standard Model. The program can generate all processes with 10 or more final state 
particles with full off-shell and finite width effects taking into account naturally both, 
spin and color correlations. The integration over the fractions xi and X2 of the initial 
partons is done via Parni [47] . 

Virtual corrections are obtained using the Helac-ILoop program [48], based on the 
Ossola-Papadopoulos-Pittau (OPP) reduction technique [49,50] and the reduction code 
CutTools [51-53]. Moreover, OneLOop [48] hbrary has been used for the evaluation of 
the scalar integrals. Reweighting techniques, and helicity and colour sampling methods 
are used in order to optimize the performance of our system. The OPP reduction at 
the integrand level takes advantage of the knowledge that the final answer for one loop 
amplitudes can be expressed in terms of a basis of known 4—, 3—, 2— and 1— point scalar 
integrals: boxes, triangles, bubbles and tadpoles(^): 

^ = ^ d,// + ^ cjf + ^ b,if + ^ ajl + n 



where TZ is the so called rational part and di,bi,Ci, ai are coefficients which have to be de- 
rived. The OPP method aims at computing them directly avoiding any computationally 
intensive integral reduction. 

The OPP reduction is based on a representation of the numerator of amplitudes, a 
polynomial in the integration momentum, in a basis of polynomials given by products 
of the functions in the denominators. Clearly, the cancellation of such terms with the 
actual denominators will lead to scalar functions with a lower number of denominators. 
By virtue of the proof provided by the Passarino-Veltman reduction [54] , we will end up 
with a tower of four-point and lower functions, as mentioned before. The determination of 
the decomposition in the new basis proceeds recursively, by setting chosen denominators 
on-shell. This is where the OPP method resembles generalized unitarity [55-64]. For 
most recent applications see 6.17. [65-70]. It is important to stress, that working around 
four dimensions, allows to compute the numerator function in four dimensions. The 
difference to the complete result is of order e, and can therefore be determined a posteriori 
in a simplified manner [52,53]. Since the calculation of the coefficients of the reduction 
requires the evaluation of the numerator function for a given value of the loop momentum, 
the corresponding diagrams can be thought of as tree level (all momenta are fixed) graphs. 
To complete the analogy, one needs to chose a propagator and consider it as cut. At 

(}) Tadpole integrals are present only when there are internal massive propagators. 
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Table I. - Integrated cross section at LO and NLO for tibb production at the LHC. The two 
NLO results refer to different values of the dipole phase space cutoff amax- The scale choice is 
fiR = ^if = mtop. 

„LO rrvi „NLO rn 1 „NLO rrni 

1489.2 ± 0.9 2642 ± 3 2636 ± 3 



this point the original amphtude for an n particle process becomes a tree level amplitude 
for an n + 2 particle process. The advantage is that its value can be obtained by a 
tree level automate such as Helac-Phegas. The bookkeeping necessary for a practical 
implementation is managed by a new software, Helac-ILoop. 

The OPP method has already been successfully applied to a large number of pro- 
cesses, apart from already mentioned ttbb, ttH — >■ ttbb and ttjj also to the production of 
three vector bosons, namely ZZZ, W+W^Z, W+ZZ and W+W~W+ final states at the 
LHC [71] and to the calculation of one-loop QED corrections to the hard-brcmsstrahlung 
process e~e+7 at e~e+ colliders [72]. Recently the OPP-approach has been implemented 
in the another framework called Samurai [73] , together with an extcntion which accom- 
modate an implementation of the generalized d-dimcnsional unitarity-cuts technique. 

The singularities from soft or coUincar parton emission arc isolated via dipole subtrac- 
tion for NLO QCD calculations [74] using the formulation for massive quarks [75] and for 
arbitrary polarizations [76]. After combining virtual and real corrections, singularities 
connected to collinear configurations in the final state as well as soft divergencies in the 
initial and final states cancel for infrared-safe observables automatically. Singularities 
connected to collinear initial-state splittings are removed via factorization by PDF redef- 
initions. Calculations are performed with the help of the Helac-Dipoles software [76], 
which is a complete and publicly available automatic implementation of Catani-Seymour 
dipole subtraction and consists of phase space integration of subtracted real radiation 
and integrated dipolcs in both massless and massive cases. The phase space restriction 
on the contribution of the dipolcs as originally proposed in [77, 78] is also implemented. 
Two values of the unphysical cutoff are always considered; amax = 1, which corresponds 
to the case when all dipolcs are included, and amax = 0.01. The independence of the 
final result on this cutoff is explicitly checked in all our results, both for the integrated 
cross section and for the differential distributions. Moreover, also in this part helicity 
sampling methods are used in order to speed up the calculation. 

The cancellation of divergences between the real and virtual corrections is always 
verified. In addition, the numerical precision of the latter has been assured by using 
gauge invariance tests and use of quadruple precision. Let us emphasise that all parts 
are calculated fully numerically in a completely automatic manner. 

Finally, the phase-space integration is performed with the multichannel Monte Carlo 
generator Phegas [41] and Kaleu [79]. 

3. Results 

We consider proton-proton collisions at the LHC with a center of mass energy of 
y/s = 14 TeV. The mass of the top quark is set to be mt = 172.6 GeV. We leave it 
on-shell with unrestricted kinematics. The jets are defined by at most two partons using 
the kr algorithm with a separation AR = 0.8, where Ai? = ^/(yi — 2/2)^ + (</>! — </'2)^; 
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Table II. - Integrated cross section at LO and NLO for ttjj production at the LHC. The two 
NLO results refer to different values of the dipole phase space cutoff amax- The scale choice is 



[Pb] 



^NLO 



:l [pb] 



^^i° =0.01 [pb] 



120.17 ± 0.08 
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Fig. 1. - Scale dependence of the total cross section for pp — !> ttbb + X (left panel) and for 
pp — ^ tijj + X (right panel) at the LHC with fiR = /xf = C ' A'o where no = mt = 172.6 GeV. 
The blue dotted curve corresponds to the LO whereas the red solid to the NLO one. 
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Fig. 2. - Distribution of the invariant mass in^j, (left panel) and the distribution in the transverse 
momentum pt^^ (right panel) of the bottom-anti-bottom pair for pp — >■ tibb + X at the LHC. 
The blue dotted curve corresponds to the LO whereas the red solid to the NLO one. 
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Fig. 3. - Distribution in the transverse momentum pt of the 1st hardest jet (left panel) and the 
2nd hardest jet (right panel) for pp — >■ tijj -I- X at the LHC. The blue dotted curve corresponds 
to the LO whereas the red solid to the NLO one. 



yi = l/2\n{Ei — pi^z) I [Ei +Pi,z) being the rapidity and 4>i the azimuthal angle of parton 
i. Moreover, the recombination is only performed if both partons satisfy |yi| < 5 (ap- 
proximate detector bounds). We further assume for tthh (tijj) processes, that the jets 
are separated by AR = 0.8 (1.0) and have |yjet| < 2.5 (4.5). Their transverse momentum 
is required to be larger than 20 (50) GeV respectively. We consistently use the CTEQ6 
set of parton distribution functions, i.e. we take CTEQ6L1 PDFs vifith a 1-loop running 
as in LO and CTEQ6M PDFs with a 2-loop running a^ at NLO. 

We begin our presentation of the final results of our analysis with a discussion of the 
total cross section. For the central value of the scale, fifi = ^p = ^q = rrit, results for 
tihh production are summarized in Tab. I whereas results for ttjj production in Tab.IL 
From the above result one can obtain K factors 



K. 



pp^ttbb+X 



= 1.77 



-''^pp-fttjj+x — 0- 



In case of pp — >■ tibb + X corrections are large of the order of 77%. However, they 
can be reduced substantially, even down to —11%, either by applying additional cuts or 
by a better choice of factorization and renormalization scales as already suggested by 
Bredenstein et al. [37]. In case oipp — >■ tijj + X we have obtained negative corrections of 
the order of 11%. In both cases a dramatic reduction of the scale uncertainty is observed 
while going from LO to NLO. The residual scale uncertainties of the NLO predictions 
for the irreducible background are at the 33% level, while for the reducible background 
the error obtained by scale variation is of the order of 11%. The scale dependence of the 
corrections for both processes is graphically presented in Fig. 1. 

While the size of the corrections to the total cross section is certainly interesting, it is 
crucial to study the corrections to the distributions. In Fig. 2 the differential distributions 
for two observablcs, namely the invariant mass and transverse momentum of the two-6-jet 
system are depicted for the pp — > ttbb + X process. Clearly, the distributions show the 
same large corrections, which turn out to be relatively constant contrary to the quark 
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induced case [34] . In Fig. 3 the transverse momentum distributions of the hardest and 
second hardest jet are shown for the pp — > ttjj + X process. Distributions demonstrate 
tiny corrections up to at least 200 GeV, which means that the size of the corrections to 
the cross section is transmitted to the distributions. On the other hand, strongly altered 
shapes are visible at high px especially in case of the first hardest jet. Let us underline 
here, that corrections to the high pT region can only be correctly described by higher 
order calculations and are not altered by soft-collinear emissions simulated by parton 
showers. 

4. Summary 

A brief summary of the calculations of NLO QCD corrections to the background 
processes pp — )- tthh + X and pp -^ ttjj + X dX the LHC has been presented. They have 
been calculated with the help of the Helac-Nlo system. 

The QCD corrections to the integrated cross section for the irreducible background 
are found to be very large, changing the LO results by about 77%. The distributions show 
the same large corrections which are relatively constant. The residual scale uncertainties 
of the NLO predictions are at the 33% level. On the other hand, the corrections to the 
reducible background with respect to LO are negative and small, reaching 11%. The 
error obtained by scale variation is of the same order. The size of the corrections to the 
cross section is transmitted to the distributions at least for the low pT region. However, 
the shapes change appreciably at high pT ■ 
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